Electromagnetic wave scattering by a superconductor 
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Abstract 

The interaction between radiation and superconductors is explored in this paper. In particular, the calculation of a plane 
standing wave scattered by an infinite cylindrical superconductor is performed by solving the Helmholtz equation in cylindrical 
coordinates. Numerical results computed up to (9(77) of Bessel functions are presented for different wavelengths showing the 
appearance of a diffraction pattern. 
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INTRODUCTION 

Since the discovery of magnetic field expulsion from 
superconductors, commonly known as Meissner effect, 
by Meissner and Ochsenfeld in 1933 [1 , the solutions 
of Maxwell equations for a superconductor in the pres- 
ence of external fields have been analyzed in numerous 
occasions for different geometries [2]-[6]. While most of 
these studies are focused on static external magnetic 
fields, the interaction between superconductors and time- 
dependent electromagnetic fields remains unexplored in 
depth in the scientific literature. 

In this paper we try to change this situation by solving 
the problem of a standing wave scattered by an infinite 
cylindrical superconductor. The problem is addressed 
by starting from the time-dependent Maxwell equation 
for the source-free case (wave equation) outside the su- 
perconducting region. For a standing wave, the space 
and time components of the wave equation solution are 
independent, reducing it to the Helmholtz and simple 
harmonic motion equations. Finally, a linear combina- 
tion of cylindrical Bessel functions (general solutions of 
the Helmholtz equation) is presented as the final solu- 
tion of the problem respecting the boundary conditions 
at infinity and at the superconductor surface. The main 
difficulties in problems like this may arise from the sym- 
metries of the system and boundary condition problem. 
Unlike the cylindrical case where the cylindrical symme- 
try simplifies the problem, in other geometries such as 
the spherical, the problem becomes non-trivial. 

To the best of our knowledge, the only extensive study 
related to the subject in the linear domain was per- 
formed by Nye in 2003 [7] where a detailed analysis of 
the scattering of plane electromagnetic waves by perfectly 
conducting regions is presented. The results are consis- 
tent with ours. Other works addressing the interaction 
between electromagnetic waves and nonlinear supercon- 
ducting materials have been presented mainly within the 
microwave community. A complete study with numerical 
computations can be found in |8j. 



STANDING WAVE SCATTERING BY AN INFI- 
NITE CYLINDRICAL SUPERCONDUCTOR 

Consider an infinite cylindrical superconductor with 
radius R in the presence of a plane standing electromag- 
netic wave and that the cylinder axis coincides with the 
z-axis (see Figure [T]). 




FIG. 1: Infinite cylinder section. 



The Maxwell equation: 
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in the source-free case and assuming the relations B = 
V X A and E -V^ - becomes: 



(2) 



which can be simplified by taking the Lorentz gauge (V • 
0): 
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the wave equation. If the electromagnetic wave propa- 
gates from infinity in the x-axis direction with the elec- 
tric field pointing in the direction of the cylinder axis, the 
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z-axis component of the vector potential A is the only rel- 
evant component needed to compute the electromagnetic 
field: 
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As the superconducting cylinder is in the presence of 
a plane standing electromagnetic wave, the component 
Az in cylindrical coordinates (since the cylinder is infi- 
nite there are no dependence on the z-coordinate) can be 
separated into: Az = F{r^ (j))T{t)^ a space-like function 
F(r, (j)) and a time-like function T{t). This separation of 
variables gives rise to two equations: 



V^F(r,(/)) + ^ 



F{r,^)=0 



and 
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the Helmholtz equation |5] and the simple harmonic mo- 
tion equation |6] Therefore, the most general solution for 
r > R is: 

Az{r, (/), t) Em=o [^rn cos(m(/)) + sin(m(/))] 
X [CmJm{kr) + DmYm{kr)] 
X [E cos{ujt) + F sm{ujt)] . (7) 

Jm{kr) and Ym{kr) are the cylindrical Bessel functions 
and /c = where A is the wavelength. The next step is 
to find the particular solution that satisfies the boundary 
conditions of a plane standing wave coming from infinity 
towards the superconducting cylinder. 

The superconducting cylinder is considered to be of 
type I with zero penetration length, a reasonable ap- 
proximation if the cylinder radius is much larger than 
the penetration length. The magnetic field is expelled 
from the superconducting region by the surface current 
defined in the coupling equation between the super- 
conductor and the magnetic field: 



-— n X B+, 
47r 



(8) 



where h is the unitary vector orthogonal to the supercon- 
ductor surface and B+ is the magnetic field in the outer 
region of the superconductor surface. As the magnetic 
field and currents are zero inside the superconducting 
cylinder, the component Az must be constant inside (set 
to zero for convenience). As there cannot be any an- 
gular dependencies on the surface of the cylinder, then 
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and: 



Az{r, 0, t) Z)m=o cos(m(/)) + Bm sin(m(/))] 

X [CraJraikr) — Cm L -r!\^ra{kr)\ 

Ym{kR) 

X [E cos(cjt) + F sm{ujt)] . (9) 



At infinity (r oo), the solution must take the form 
of the plane standing wave assumed to be: 



Az{r, ^, t) = Aq cos{kr cos (j)) cos(a;t). 



(10) 



where Aq is the amphtude of the wave and x = rcos( 
According to Jacobi- Anger expansion 



cos(2; cos <p) = Jo{z) + 2 (— l)"'^2m(-z) cos(2m0), 



and therefore the final solution reads for r > R: 
F{r,(P] 
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and F(r, ^) = for r < R. Such a mathematical result is 
hard to interpret and visualize but possible to represent 
with enough accuracy through numerical computations 
up to m = 77. The spatial representation and analy- 
sis of Az for different wavelengths are presented in the 
following section. 



SIMULATION 

The numerical computations that simulate the scatter- 
ing were performed in a 2D array representing the cross- 
section of the superconducting cylinder (centered at the 
origin) in the x-y plane ranging from minus 40 to 40 in 
steps of 0.1 radius units for both x and y directions at 
the instant t = 0. The simulation of the plane stand- 
ing wave scattering by the superconducting cylinder is 
shown in Figure [2] in arbitrary units for A = AR (top) and 
A = 6R (bottom). For A = AR the scattering presents 
a diffraction pattern of nodes and anti-nodes around the 
superconducting cylinder while converging to the plane 
wave at large distances. The scattering effects are per- 
haps more visible for A = 6i?, the nodes and anti-nodes 
around the cylinder are clearly amplified and separated 
by roughly the wavelength. The nodes and anti-nodes are 
maximal near the superconductor and decrease gradually 
with distance. 



DISCUSSION 

According to Maxwell equations, the interaction be- 
tween a plane wave and a superconductor gives rise to 
a diffraction pattern. We would like to challenge exper- 
imental groups to test the result obtained in this paper 
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even though we are aware of some technical hmitations 
that may constrain the observations. For example, a dif- 
ficulty may come from generating a large enough con- 
trolled plane standing wave with a wavelength of the 
same order of magnitude as the superconductor radius 
to enhance the scattering. Also, it is important to stress 
that the wave period must be much larger than the super- 
conductor relaxation time so that it can respond instantly 
to field variations. 
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FIG. 2: Standing wave scattered around the infinite cylindri- 
cal superconductor for A = 4i? (top) and A = 6R (bottom). 
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